Preliminary heavy-light decay constants from the MILC collaboration* 
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Preliminary results from the MILC collaboration for fs, /s s , /z>, /d„ and their ratios are presented. We 
compute in the quenched approximation at fi = 6.3, 6.0 and 5.7 with Wilson light quarks and static and Wilson 
heavy quarks. We attempt to quantify all systematic errors other than quenching. 



1. PRELIMINARIES 

Over the past year, we have been computing 
heavy-light decay constants in the quenched ap- 
proximation on Intel Paragon computers. Most 
of the computations have been performed on the 
512-node Paragon at Oak Ridge National Labora- 
tory, but Paragons at Indiana University and at 
the San Diego Supercomputer Center have also 
been used. In many respects the calculations are 
standard, and we emphasize here only the distin- 
guishing features. 

The initially very slow I/O speeds of the 
Paragon and the lack of long-term storage capa- 
bility at ORNL forced us to do all the computa- 
tions "on the fly." The hopping parameter com- 
putation of the heavy quark propagator, as sug- 
gested by Henty and Kenway [1], makes such on- 
the-fly computations possible for heavy-light sys- 
tems. Because the full light and heavy propaga- 



Table 1 

Lattice parameters. 



name 


P 


size 


# 


configs. (planned) 


A 


5.7 


8 3 x 48 




200 (200) 


B 


5.7 


16 3 x 48 




100 (100) 


C 


6.0 


16 3 x 48 




48 (100) 


D 


6.3 


24 3 x 80 




98 (100) 



* presented by C. Bernard 



tors for all spin-color sources can not be stored in 
memory, we work only with one spin-color source 
for light and heavy at a time, and restrict our- 
selves to mesons which are diagonal in spin-color 
(i.e., pseudoscalars and the z-component of vec- 
tors). We run 400 hopping parameter passes. At 
f3 = 6.3, this gives very good convergence for 
heavy quarks with n^ < 0.145. 

Gaussian quark sources in Coulomb gauge are 
used. The overrelaxed gauge fixer is run until the 
sum of the trace of all spacelike links (normalized 
to 1 when all links are unit matrices) changes by 
less than 7 x 10" 7 per pass. On the 24 3 x 80 



lattices at /? = 6.3, this takes about 435 passes. 
The half-width of the gaussian is ps 0.4 fm. 

We compute "smeared-local" and "smeared- 
smeared" propagators. Because the mesons must 
be constructed at each of the 400 orders of the 
hopping parameter expansion, it is too expen- 
sive to sum the central point of the smeared 
sinks over the entire spatial volume, even using 
FFT's. Instead, we simply sum over a subset 
of the points in the spatial volume. This allows 
intermediate states of non-zero 3-momentum to 
contribute. For the heavy-light mesons studied 
here, the higher momentum states are well sup- 
pressed at asymptotic time by their higher en- 
ergy. However, the static-light mesons have no 
such suppression, and the contribution of higher 
momentum states is limited only by their overlap 
with the sources. 

We sum the sinks over 16 points on a time slice. 
Using computed static-light wavefunctions [2], we 
find that the contamination in static-light decay 
constants from nonzero momentum states is small 
(ps 0.7%) for lattices with spatial size of ps 1.5 fm 
(lattices A, C, and D). However, on lattice B, 
with spatial size of ps 3 fm, the contamination is 
p» 60%. We thus do not include the static point 
from lattice B (nor from lattice A, so we may 
compare the /? = 5.7 lattices without bias). 

Since we only have results for degenerate light 
quarks, we determine k s , the strange quark hop- 
ping parameter, by adjusting the pseudoscalar 
mass to \j2m 2 K — m\, the lowest order chiral per- 
turbation theory value. 

For heavy-light mesons we use the Kronfeld- 
Mackenzie [3] norm (yl — 6k) and adjust the 
measured meson pole mass upward by the differ- 
ence between the heavy quark pole mass ("mi") 
and the heavy quark dynamical mass ("012") as 
calculated in the tadpole-improved tree approxi- 
mation [3]. 

2. RESULTS 

A plot of fp\/Mp vs. 1/Mp is shown in fig. 1 
for lattice D. The fit is covariant, to the form 
c + ci/Mp + c 2 /Mp. The x 2 / d -°- f for the nt is 
ps 2 (confidence level ps 10%), whether or not the 
static-light point is included. The rather low con- 
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Figure 1. f P (M P )- vs. 1/M P for lattice D. The 
light quark is extrapolated to the physical mass 
(m„ +m d )/2. 



Table 2 

Results for decay constants 

MeV scale used throughout. 



and ratios. / 7r = 131 





A 


B 


C 


D 


!b 


195(6) 


198(4) 


175(5) 


166(4) 


Sb. 


244(5) 


237(3) 


207(4) 


192(3) 


Id 


227(5) 


227(4) 


205(4) 


198(2) 


Id, 


275(4) 


273(3) 


239(3) 


225(2) 


Sb, 
!b 


1.25(2) 


1.20(1) 


1.18(1) 


1.16(1) 


/d, 
Id 


1.21(1) 


1.20(1) 


1.17(1) 


1.13(1) 



fidence level may perhaps be due to the fact that 
we have not included additional large-ma correc- 
tions to the action and operators [4], or simply 
to the small differences between the heavy quark 
mass and the meson mass Mp. Such effects are 
under investigation. Note that, in an earlier cal- 
culation [5], the statistical errors were consider- 
ably larger, and the x 2 /d.o.f for such fits was 
good. Here the level of statistical precision has 
increased to a level where small effects are be- 
coming important. 

Table 2 shows results from the four lattices. 
The lattice-spacing dependence is apparent, but 
little, if any, finite volume effect is present (com- 
pare A and B). This is seen more clearly in Fig. 2, 
which shows fp vs. lattice spacing. It is natu- 
ral to extrapolate the /^-scale results linearly to 
the continuum; we get 147(6) MeV. Note that 
the /^-scale results have much less a dependence 
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Figure 2. fg vs. a. Diamonds have scale set 
by f w ; crosses, by m p . Fit is to the diamonds; 
"fancy diamond" gives the a = extrapolation. 
The higher cross at a ps 0.7 is from lattice B. 

than those using m p . This makes sense since f w 
and fs are likely to have rather similar finite a 
effects. If the m^-scale results are extrapolated 
linearly to a = 0, the result is 119(8) MeV, con- 
siderably smaller than the /^-scale extrapolation. 
However, there also seem to be larger finite vol- 
ume effects in the m^-scale results, which is rea- 
sonable since the p is a larger state than the it. 
If we first adjust upward the results from lattices 
C and D by the presumed finite volume correc- 
tion obtained by comparing the lattices B and A, 
the result of the m^-scale extrapolation (144(9) 
MeV) is consistent with the /^-scale result. 

The results on lattice D are consistent with 
those of [5]. The major cause of the somewhat 
smaller central values here is a lower (but still 
consistent) estimate of the scale (1/a ps 3.0 GeV 
here vs. ps 3.2 GeV in [5]). 

We linearly extrapolate to a = all results in 
Table 2. Systematic errors are then estimated — 
in a very preliminary fashion — as follows: (1) 
Changes of fitting ranges (in t) for the propaga- 
tors and of types of fits in \/M for fp\/Mp give 
a typical variation of about twice the statistical 
errors. (2) The dependence on the determination 
of k s is estimated by finding the change in the 
extrapolated results if k s is fixed instead using 
the vector state <f>. The difference is especially 
significant for JbJJb and is ps 0.05 there. (3) 
Finite volume effects are estimated by taking the 



fractional difference between results from lattices 
A and B, using the m p scale. This is conserva- 
tive, since the central values are found with the 
f w scale. (4) We estimate scale errors by compar- 
ing the results at /? = 6.3 with f w and m p scales. 
The difference (ps 13 MeV for the decay constants 
and ps 0.02 for the ratios) is roughly comparable 
to what we would get by comparing extrapolated 
f w - and nip-scale values, after adjusting for the 
apparent finite volume effects. (5) The effects 
of using heavy Wilson fermions without the ad- 
ditional corrections to the action and operators 
detailed in [3,4] are estimated by comparing the 
original fits at /? = 6.3 (see, e.g., Fig. 1) with 
fits using only the 6 lightest heavy-light states 
(and the static-light point). In the original fits 
the maximum value of {m^ — m\)jm2 is 0.22; in 
the new ones, 0.04. The differences in the results 
are quite small: ps 4 MeV for the decay constants 
and ps 0.01 for the ratios. 

Adding the above systematic errors in quadra- 
ture, our preliminary results are 

f B = 147(6)(23); /d = 181(4)(18); ; (1) 

f B . = 164(5)(20); f D . = 195(3)(16); (2) 



- 1-13(2)(8); 
JB Id 



1.09(1)(4), (3) 



where the decay constants are in MeV. Study of 
the quenching errors is in progress. 
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